第 36 卷 第 11 期 计算 机 应 用 研究 36No 
录用 定稿 Application Research of Computers Accepted Paper 


光谱 信息 ， 拥 有 数量 庞大 的 细微 光谱 特征 ， 因 此 解决 了 光谱 信 ” 据 的 稀 玻 表示 系数 构造 相似 度 和 矩阵 ， 然 后 利用 谱 聚 类 方法 得 到 
息 不 足 的 问题 ， 提 高 了 地 物 目标 的 区 分 和 识别 能 力 由 。 但 由 于 数据 的 子 空间 聚 类 结果 ， 然 而 稀 玻 子 空间 聚 类 在 聚 类 过 程 中 却 


相对 于 传统 的 多 光谱 数据 ， 高 光谱 遥感 影像 提供 了 丰 


融合 高 光谱 影像 三 维 空 谱 特征 的 子 空间 聚 类 算法 


李 文 洲 *， 邓 秀 勤 *， 刘 语 春 ， 
(广东 工业 大 学 a. 应 用 数学 学 院 ; b. 计算 机 学 院 ， 广州 510006) 


摘 要 : 为 提高 高 光谱 喀 感 影像 的 聚 类 精度 ， 将 三 维 空 谱 特征 和 子 空间 聚 类 算法 相 结 合 ， 提 出 一 种 新 的 稀 政 子 空间 聚 
类 模型 ， 在 关注 高 光谱 影像 的 光谱 信息 的 同时 ， 也 关注 空间 上 下 文 信息 。 首 先 提取 高 光谱 影像 像素 点 的 三 种 三 维 空 谱 
特征 ， 然 后 通过 特征 对 子 空间 聚 类 模型 的 系数 算 阵 进行 加 权 ， 使 得 像素 点 可 被 与 它 最 为 相似 的 像素 点 稀 足 表示 ， 从 而 
获得 更 好 的 系数 答 阵 ， 最 后 由 系数 和 矩 阵 通 过 谱 聚 类 获得 更 好 的 聚 类 结果 。 算 法 对 四 个 经 典 高 光谱 数据 集 进行 实验 ， 并 
将 实验 结果 与 六 种 聚 类 算法 进行 比较 , 结果 表明 ,所 提出 的 3DF-SSC 算法 在 四 个 数据 集 上 获得 的 聚 类 精度 都 比 其 他 算 
法 要 高 ， 对 于 同样 是 利用 三 维 空 谱 特征 的 M3DF3 算法 ，3DF-SSC 算法 最 高 能 提高 8.62% 的 精度 ， 而 与 同样 是 利用 空间 
上 下 文 信息 对 子 空 间 聚 类 算法 进行 改进 的 L2-SSC 算法 和 SS-LRSC 算法 相 比 ， 最 高 能 提高 25.18% 的 精度 。 
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Subspace clustering algorithm fused three-dimensional 
spatial spectral features of hyperspectral Images 


Li Wenzhou’, Deng Xiuqinsi, Liu Fuchuny 
(a. School of Applied Mathematics, b. School of Computers Guangdong University of Technology, Guangzhou 510006, China) 


Abstract: In order to improve the clustering accuracy of hyperspectral images, this paper proposed a new sparse subspace 
clustering model combined three-dimensional spatial spectral features with subspace clustering algorithms. While focusing on 
the spectral information of hyperspectral images, it also paid attention to spatial context information. Firstly, extracting three 
kinds of three-dimensional spatial spectral features from the pixels of the hyperspectral image. Then the features weighted the 
coefficient matrix of the subspace clustering model so that the pixel points could sparsely represent the pixel point to which they 
are most similar, thereby obtaining the better coefficient matrix. Finally, it used the coefficient matrix to obtain better clustering 
results with spectral clustering. The algorithm experimented on four classical hyperspectral datasets, and compared the 
experimental results with six clustering algorithms. The results show that the proposed algorithm achieves higher clustering 
accuracy on the four datasets than the other algorithms. The algorithm can achieve at most 8.62% accuracy than the algorithms 
based 3D spatial spectral features like M3DF’ algotithm and 3DF-SSC algorithm, and at most 25.18% than the algorithms which 
improve the subspace clustering algorithm by using spatial context information like L2-SSC algorithm and SS-LRSC algorithm. 


Key words: hyperspectral images; three-dimensional spatial spectral features; subspace clustering; spatial context information 


近年 来 , 稀 疏 子 空间 限 类 (sparse subspace clustering, SSC) 算 


引言 


法 外 也 被 应 用 于 高 光谱 遥感 影像 的 聚 类 中 ， 并 显示 


上 上 


巨大 的 洪 


的 力 。SSC 算法 是 一 种 基于 谱 聚 类 的 数据 聚 类 框架 ， 它 可 利 


在 监督 分 类 中 ， 对 高 维 高 光谱 遥感 影像 的 大 量 训练 样本 进行 标 ” 仅 利 用 了 光谱 信息 ， 没 有 考虑 空间 上 下 文 信息 ， 导 致 其 最 终 的 
记 具 有 一 定 的 困难 ， 因 而 K-means 聚 类 请 和 谐 聚 类 (SC)PI 等 无 到 类 图 像 会 出 现 大 量 椒盐 噪声 钻 。 对 此 ， Zhai 等 人 提出 L2- 
监督 聚 类 算法 被 广泛 应 用 于 高 光谱 遥感 影像 的 聚 类 中 SSC 算法 , 融入 邻 域 信 息 对 系数 进行 TV (total variation ) 约束 ， 
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从 而 提高 聚 类 精度 。 为 了 强化 高 光谱 遥感 影像 的 空间 信息 , Zhu 。” 其 中 多 和 疙 是 两 个 数据 样本 ， 建 立 加权 稀 疏 表示 模型 ; 
等 们 提出 高 光谱 影像 的 三 维 空 谱 特 征 提取 方法 ， 并 基于 稀 玻 3 1 ca 
min —IC;, 
决策 融合 算法 提出 M3DF 算法 对 高 光谱 影像 进行 地 物 分 类 ,而 EAG (4) 


Xu 等 人 辐 提 


出 的 SS-LRSC 算法 则 是 利用 空间 信息 对 子 空间 聚 


类 算法 求 得 的 系数 矩阵 进行 调制 , 从 而 获得 了 


为 了 更 有 效 地 利 
文 提出 一 种 融合 高 光谱 影 


j 高 光谱 遥感 影像 的 空间 上 下 


更 高 的 聚 类 精度 。 


像 三 维 空 谱 特 征 的 子 空 


文 信息 ， 本 
间 聚 类 算法 


(three dimension spatial feature-Sparse Subspace clustering, 3DF- 


SSC)， 通 过 提取 像素 的 三 维 空 谱 特 征求 出 像素 之 间 的 相似 度 ， 


相似 的 像素 


尽 可 
性 。 实验 结 果 表 
高 光谱 数据 身 


能 参与 到 数据 的 自 表示 当中 ， 


都 能 获得 准确 可 靠 的 聚 类 结果 。 


A 


1 ” 稀 琉 子 空间 聚 类 


黎 琉 子 


2 
i 
基本 思想 是 对 给 定 的 


间 聚 类 是 一 种 基 


明 ， 本文 提 出 的 3DF-SSC 算法 对 于 四 


于 稀疏 表示 和 谱 聚 


在 子 空间 聚 类 的 稀 玻 度量 中 引入 与 相似 度 成 反比 的 权重 ， 迫 使 


提高 聚 类 的 准 
个 经 


类 方法 ， 


系数 ， 由 此 构造 相 
聚 类 结果 。 


高 光谱 遥感 


中 得 到 数 


全 


中 ， 首 先 对 像素 按 波段 D 进行 列 化 ， 将 三 


组 高 维 数据 ， 


以 度 和 矩阵 ， 再 利 | 


确 
的 


类 的 子 空间 聚 
寻找 其 稀 琉 表示 
j 谱 聚 类 方法 ， 如 规范 化 割 


影像 (hyperspectral image, HSI) J eR”N*? 


射 到 二 维 


空间 了 [yy yyw] ER 9 然后 假设 来 


有 像素 都 是 属 了 


的 高 光谱 影像 了 映 


了 的 所 


F 1 个子 空间 的 并 Ui15; ,每 个 子 空间 对 应 一 个 特 


定 的 陆地 覆盖 类 型 ， 最 后 子 空 间 聚 类 利 


将 每 个 像素 表示 为 其 他 像素 的 线 怕 


为 不 同 的 类 。 


其 中 : Ye R?* 为 二 维 囊 S7 入 


E Ee€ RP 为 


multipliers) 


在 MATLAB 
法 求 


具体 模型 为 


minlCl + 4lel; 


St. 


误差 矩阵 ，4 > 0 为 平衡 参数 。 


组 合 ， 从 而 


解 模型 式 (1) 可 得 稀 玻 表示 系数 
矩阵 C 通过 式 (2) 构 造 邻接 矩阵 We RY ，W 
示 两 个 像素 之 间 的 相似 度 [n021。 


而 将 这 些 


Y=YC+E,diag(C)=0 


邻接 矩 隐 


自我 表示 
像素 聚 类 


(D) 


E 阵 ，M、N、DD 分别 为 高 光谱 影 
像 的 长 、 宽 和 波段 数 ，C e RW*W 为 稀 朴 表示 系数 珑 阵 ， 


使 用 ADMM(alternating direction method of 
E 阵 C ， 然 后 
的 元 素 表 


O) 
EW 得 到 


三 =|Cl+|c| 
最 后 通过 Ncut(normalized cuts) 算 法 [13] 
最 终 的 聚 类 结果 


加 权 稀 疏 子 空间 表示 模型 ， 引 入 权重 


使 得 数 和 


间 聚 类 模型 对 噪声 较为 敏感 ， 所 以 文献 [14] 提 出 


居 尽 可 能 


与 它 最 


为 相似 的 数 


加 权 ， 即 式 (3)。 


居 线 性 表示 ， 首 先 利用 


A 


2 
oo | 


斯 相似 度 函数 对 系数 进行 


SI. 


Y=YC+E,diag(C)=0 


2 高 光谱 遥感 影像 的 三 维 空 谱 特征 


高 光谱 遥感 影像 包含 了 关于 地 表 物 质 的 丰富 空间 、 辐 射 和 
光谱 三 重信 息 ， 因 此 近年 来 也 有 许多 高 光谱 遥感 影像 的 特征 提 
取 方 法 ， 但 是 现 有 的 空 谱 特征 提取 方法 大 多 是 将 不 同 波段 的 空 
间 特 征 机 械 地 组 合 起 来 ， 为 了 充分 利用 高 光谱 遥感 影像 空 谱 结 
构 中 的 上 下 文 信息 ， 文 献 [15] 提 出 了 高 光谱 遥感 影像 的 三 维 空 
谱 特 征 的 提取 方法 ， 如 三 维 LBP 空 谱 特 征 、 三 维 空 谱 开 特征 和 
维 空 谱 闭 特征 。 具 体 特征 提取 公式 如 下 : 

给 定 高 光谱 遥感 影像 HSI: Te RW™*? ,像素 点 k 的 3DLBP 


空 谱 特 征 为 


3DLBP - 名 —1) if T(3DLBP,)<P 
P+l 其 它 
其 中 

= 1 lxl20o 
:= | 对 < (0) 
T(3DLBP,) = pa 1)—s(1,—1) (7) 
其 中 : Z(i=0,…,P--D) 是 以 像素 点 为 中 心 的 外 接 球 半 径 为 1 
的 正八 面体 的 第 i 个 采样 点 的 像素 值 ，P 为 采样 点 的 数量 。 为 
了 实现 灰 度 不 变性 ， 通 过 八 面体 邻 域 的 中 心 点 的 像素 值 对 采样 


点 执行 二 值 化 ， 从 而 获取 局 部 空 谱 结构 的 三 维 局 部 二 值 模式 ， 
如 式 (5) 所 示 ， 其 中 元 为 中 心 点 ，Ga 为 判断 冰 值 。 即 若 顶 点 位 
的 像素 值 与 中 心 点 的 像素 值 的 差 xX 大 于 o ， 将 其 设 为 1， 否 贝 
设 为 0。 


i 


像素 点 k 的 三 维 空 谱 开 特征 为 

3DOR 二 max(Gmin(7(2 》, z) E CGOm， yo; Zz0))) (8) 
像素 点 上 的 三 维 空 谱 闭 特征 为 

3DCR = min(max(1(%, 了》， z) E O(N, y0， Zz0))) (9) 


其 中 : 7(%w, yozo) 是 位 于 (xm,yozo) 处 的 像素 值 。Q(%o,y0,z0) 是 
高 光谱 影像 HSI 中 以 像素 点 (xm, yzo) 为 中 心 的 立方 体 集合 。 

3 ”融合 高 光谱 三 维 空 谱 特征 的 子 空间 聚 类 算法 

由 于 基于 子 空间 聚 类 算法 的 高 光谱 遥感 影像 的 地 物 区 分 仅 
了 影像 的 光谱 信息 ， 因 此 本 文 将 三 维 空 谱 特 征 与 子 空间 聚 
类 算法 进行 融合 ， 利 用 影像 的 光谱 信息 的 同时 ， 还 利用 像素 的 
空间 上 下 文 信息 ， 使 得 每 个 像素 都 可 以 被 与 它 最 为 相似 的 像素 
稀 玻 表示 ， 提 高 子 空 间 聚 类 精度 。 


首先 对 高 光谱 遥感 影像 中 的 每 个 像素 点 ,量化 以 


利用 


像素 点 
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k 为 中 心 的 三 维 局 部 正八 面体 领域 的 空 谱 结构 ， 通 过 八 面体 邻 


前 分 别 取 4 为 0.001、0.005、0.01、0.05、0.1 或 0.5， 对 数据 集 


域 的 中 心 点 的 像素 值 对 采样 点 执行 二 值 化 ， 从 而 统计 每 个 像素 
点 大 的 八 面体 邻 域内 的 3DLBP 特征 : 3DLA.， 然后 对 以 像素 点 
k 为 中 心 的 立方 体 集合 进行 膨胀 、 腐 蚀 、 开 启 和 闭合 运算 ， 提 
取 每 个 像素 点 的 三 维 空 谱 开 特 征 和 闭 特征 : 3DOR 、3DCR 。 


寻找 对 模型 较 优 的 参数 。 实 验 结果 表明 当 4 取 0.01 时 模型 能 得 
到 较 高 的 聚 类 精度 。 

将 本 文 提出 的 3DF-SSC 算法 的 实验 结果 与 -meansD、 
SCBI]、SSCm、M3DF3I、SS-LRSC 算法 [I 和 L2-SSCIG 算 法 的 


将 三 种 特征 分 别 记 做 X\、，X>、X3， 通过 高 斯 函数 计算 每 两 个 像 

素 点 之 间 的 相似 度 为 : 

COC ed (C.F. He. Pe. 
re 


Gy = exp( | 


(10) 
其 中 : (X,), 和 (X,), 分 别 为 像素 了 和 9 的 第 n 个 特征 ， 
n=1,2,3。 
最 后 用 1/G; 对 系数 的 4 范 数 进行 加 权 , 使 得 相似 度 越 大 的 
像素 点 权重 越 小 ， 相 似 度 越 小 的 像素 点 权重 越 大 ， 让 像素 点 尽 
可 能 被 最 相似 的 像素 线性 表示 ， 具 体 模型 如 下 : 


1 
min > 一 -|Cil+4|ZcC 一 2 
Sa (11) 
st. C=0,1=1,2,..,MN 
3DF-SSC 算法 步骤 如 下 : 
输入 : HSI:TeR”N?， 聚 类 数 1。 
a) 对 高 光谱 影像 了 进行 降 维 ， 得 Ye R?*; 


b) 提 取 Y 中 每 个 像素 的 三 种 特征 , 3DLBP: X =[Xx4,…,Xip]、 


3DOP: X， =[X%1,*…*, Xp] 和 3DCP: X3 三 [xp] 六 


co) 由 式 (10) 计 算 像 素 之 间 的 相似 度 ; 
dj) 求解 模型 (11), 得 子 空间 表示 系数 矩阵 C ; 
6) 由 式 (2) 计 算 相 似 度 矩阵 W ,并 利用 Neut 算法 得 到 聚 类 结 


输出 : 图 像 聚 类 结果 。 
4 。 实验 结果 与 分 析 
为 了 评估 算法 的 聚 类 精度 ， 本 文 将 对 四 个 经 典 高 光谱 数据 


集 : Salinas-A、the Pavia Center、 Pavia University 和 Indian Pines09] 
进行 实验 测试 。Salinas-A 高 光谱 遥感 影像 由 AVIRIS 传感器 在 
California 的 Salinas Valley 获取 , 数据 集 包 含 224 个 波段 图 像 ， 
每 个 波段 图 像 大 小 为 86*83， 分 为 6 类 地 物 ; Pavia Center 高 光 
谱 遥 感 影像 由 ROSIS-03 传感器 在 意大利 的 由 维 亚 中 心 获取 ， 
简称 PaviaC。 数 据 集 包 含 102 个 波段 图 像 ， 每 个 波段 图 像 大 小 
为 1096*715， 总 共 148152 个 标记 样本 ,分 为 9 类 地 物 。 实 验 
采用 其 中 尺寸 为 120*120， 包 含 6 类 地 物 的 高 光谱 图 像 进行 测 
试 ; Pavia University 高 光谱 遥感 影像 由 ROSIS-3 传感器 采集 于 
意大利 帕 维 亚 大 学 周边 的 帕 维 亚 大 学 ， 拥 有 103 个 波段 图 像 ， 
每 个 波段 图 像 大 小 为 610*340， 分 为 9 类 地 物 。 实 验 采 用 其 中 
尺寸 为 200*100, 包含 8 类 地 物 的 高 光谱 图 像 进 行 测试 Indian 
Pines 是 AVIRIS 传感器 于 印第安 纳 洲 的 农场 采集 ， 其 中 包含 
220 个 尺寸 为 145*145 的 高 光谱 图 像 ， 分 为 6 类 地 物 。 

由 于 模型 式 (11) 中 存在 平衡 参数 4 ， 因 此 在 进行 对 比 实验 


一 | 


实验 结果 进行 对 比 。 七 种 算法 对 四 个 数据 集 的 聚 类 精度 如 表 1 
所 示 , 其 中 PaviaU 和 SalinasA 数据 集 的 聚 类 结果 如 图 1 和 图 2 
所 示 。 K-means 和 SC 算法 是 经 典 的 无 监督 聚 类 算法 ， 从 表 1 
以 及 图 1、2 可 以 看 出 , 这 两 种 仅 关 注 光谱 信息 的 算法 取得 的 聚 
类 效果 较为 一 般 ， 对 高 光谱 像素 的 判别 能 力 较 差 ， 存 在 大 量 的 
误 分 类 和 椒盐 噪声 。 而 由 算法 的 运行 结果 可 看 出 , 在 提取 高 光 
谱 影 像 的 三 维 空 谱 特征 后 ， 影 像 的 聚 类 效果 得 到 了 显著 提高 ， 

可 见 高 光谱 遥感 影像 的 空间 上 下 文 信息 有 助 于 对 像素 点 的 识别 ; 
对 于 子 空间 聚 类 算法 , 传统 的 SSC 算法 的 运行 结果 和 k-means、 
SC 算法 一 样 较 为 一 般 ， 而 L2-SSC 算法 和 SS-LRSC 算法 添加 
空间 上 下 文 信息 后 ， 遥 感 影像 的 聚 类 效果 则 有 了 明显 的 提高 ， 

由 此 可 见 ,获得 空间 上 下 文 信息 的 子 空 间 聚 类 的 聚 类 效果 优异 。 

表 1 七 种 算法 在 四 个 数据 集 上 的 聚 类 精度 


Table 1 Clustering accuracy of seven algorithms on four datasets 


K-means SC SSC M3DF” SS-LRSC L2-SSC 3DF-SSC 


SalinasA 65.12 69.27 66.49 81.78 75.76 75.14 90.40 
PaviaC 64.37 64.98 58.76 85.66 86.67 82.16 89.47 
PaviaU 56.69 56.68 54.39 70.44 68.99 63.30 74.15 
IndianP 34.69 34.69 37.23 68.22 46.13 44.54 69.72 


K-means SC SSC 


级 
pS 
f 


SS-LRSC 
图 1 七 种 算法 对 PaviaU 数据 集 的 聚 类 结果 


Fig.l Clustering maps of seven algorithms with PaviaU dataset 
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MI3DF” 3DF-SSC 
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后 ， 


原 图 K-means SC SSC 


M3DF’ SS-LRSC L2-SSC 3DF-SSC 
图 2 七 种 算法 对 SalinasA 数据 集 的 聚 类 结果 
Fig.2 Clustering maps of seven algorithms with SalinasA dataset 

表 1 以 及 图 1、2 的 实验 结果 表明 ， 在 融合 三 维 空 谱 特征 
本 文 提 出 的 3DF-SSC 算法 的 聚 类 效果 都 比 其 他 算法 要 好 ， 


四 个 数据 集 的 聚 类 精度 都 能 达到 70% 左 右 或 以 上 ， 尤 其 是 
SalinasA 的 聚 类 精度 高 达 90.40%。 对 于 同样 是 利用 三 维 空 谱 特 
征 的 M3DFs 算 法 ， 本 文 算法 获得 的 聚 类 精度 也 能 够 有 所 提高 ， 
最 高 能 提高 8.62%; 而 与 同样 是 利用 空间 上 下 文 信息 对 子 空间 
聚 类 算法 进行 改进 的 L2-SSC 算法 和 SS-LRSC 算法 相 比 ， 本 文 
算法 的 运行 结果 也 更 为 优异 ， 尤 其 是 IndianP 数据 集 ， 聚 类 精 


度 E 


EL2-SSC 算法 提高 了 25.18%。 由 此 可 见 ， 本文 提 出 的 3DF- 


SSC 算法 能 有 效 对 高 光谱 遥感 影像 进行 聚 类 。 
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结束 语 


| 


本 文 提出 的 结合 三 维 空 谱 特征 和 子 空 


聚 类 算法 的 3DF- 


SSC 算法 ， 在 关注 高 光谱 影像 的 光谱 信息 的 同时 ， 也 关注 空间 


上 下 文 信息 。 通 过 提取 高 光谱 影像 像素 点 的 三 种 三 维 空 谱 特征 


j 高 斯 函数 由 多 特征 对 子 空 间 聚 类 算法 的 系数 矩阵 进行 加 权 ， 


以 


建立 了 融合 三 维 空 谱 特 征 的 稀 玻 子 空间 聚 类 模型 。 模 型 对 四 种 


经 典 高 光谱 数据 入 


进行 实验 ， 并 将 实验 结果 与 六 种 经 典 隧 类 算 
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法 相 比 较 。 实 验 结果 表明 ， 本 文 提出 的 3DF-SSC 算法 在 四 个 数 
据 集 上 都 获得 了 良好 的 聚 类 精度 ， 并 且 比 其 他 经 典 算 法 的 实验 


一 


结果 都 要 更 好 。 
参考 文献 : 


[1] 


[2] 


[3] 


邓 贤 明 , 苗 放 ， 翟 涌 光 , 等 . 基于 形态 学 的 两 种 高 光谱 目标 探测 改进 算 
法 [J]. 中 山大 学 学 报 : 自然 科学 版 , 2017, 56 (1): 151-160. (Deng 
Xianming, Miao Fang, Zhai Yongguang, et al. Two modified target detection 
algorithms based on morphology for hyperspectral imagery [J]. Acta 
Scientiarum Naturalum Universitatis Sunyatseni, 2017, 56 (1): 151-160. ) 
Lloyd S. Least squares quantization in PCM [J]. IEEE Trans on Information 
Theory, 1982, 28 (2): 129-137. 


Luxburg U. A tutorial on spectral clustering [J]. Statistics and Computing. 


i | ChinaXiv 合 作 期 刊 
维 空 谱 特 征 的 子 空 间 聚 类 算法 第 36 卷 


2007, 17 (4): 395-416. 

[4] Ehsan E, Ren V. Sparse Subspace clustering [C]/ IEEE Conference on 
Computer Vision and Pattern Recognition. Washington DC: IEEE Computer 
Society, 2009: 2790-2797. 

[5] 王 卫 卫 , 李小平 , 冯 象 初 , 等. 稀疏 子 空间 聚 类 综述 [可 自动 化 学 报 ， 
2015, 41 (3): 1373-1383. (Wang Weiwei, Li Xiaoping, Feng Xiangchu, et al. 
A Survey on Sparse Subspace Clustering [J]. Acta Automatica Sinica, 2015, 
41 (3): 1373-1383. ) 

[6] Zhai Han, Zhang Hongyan, Zhang Liangpei, et al. A new sparse subspace 
clusteing algorithm for hyperspectral remote sensing imagery [可 , IEEE 
Geoscience and Remote Sensing Letters, 2017, 14 (1): 43 — 47. 

[7] Zhu Jiasong, Hu Jie, Jia Sen, et al. Multiple 3D feature fusion framework 
for hyperspectral image classification [J]. IEEE Trans on Geoscience and 
Remote Sensing, 2018, 56 (4): 1873-1884. 

[8] Xu Jinhuan, Huang Nan, Xiao Liang. Spectral-spatial subspace clustering 
for hyperspectral images via modulated low-rank representation [C1]// Proc 
of IEEE International Geoscience and Remote Sensing Symposium. 2017: 
3202-3205 


[ 


bs 


Shi J, Malik J. Normalized cuts and image segmentation [J]. IEEE Trans on 

PaUern Analys is and Machine Intelligence, 2000, 22 (8): 888-905. 

[10] Manya V A, Jos M B, Figueiredo Mario A T. An augmented Lagrangian 
approach to the constrained optimization formulation of imaging inverse 
problems [J]]. IEEE Trans on Image Processing, 2011, 20 (3): 681-695. 

[11] Zhang Hongyan, Zhai Han, Zhang Liangpei, et al. Spectral-spatial 
sparsesubspace clustering for hyperspectral remote sensing images [J]. IEEE 
Trans on Geoscience Remote Sensing, 2016, 54 (6): 3672--3684. 

[12] Ehsan E, Ren V. Sparse subspace clustering: Algorithm, theory, and 
applications [J]. IEEE Trans on Pattern Analysis Machine Intelligence, 2013, 
35 (11): 2765-2781. 

[13] Yuri B, Gareth F. Graph cuts and efficient N-D image segmentation [J]. 
Computer Visual, 2006, 70 (2): 109—131. 

[14] 李涛 , 王 卫 卫 ,有 翟 栋 等 . 图 像 分 割 的 加 权 稀 疏 子 空间 聚 类 方法 [J]. 系 
统 工程 与 电子 技术 , 2014, 36 (3): 580-585. Li Tao, Wang Weiwei, Zhai 
Dong, et al. Weighted-sparse subspace clustering method for image 
segmentation [J]. Systems Engineering and Electronics, 2014, 36 (3): 580- 
585. 

[15] 胡 杰 . 高 光谱 遵 感 影像 三 维 空 谱 特征 提取 与 小 样本 分 类 技术 研究 [D]. 
深圳 : 深圳 大 学 , 2017. (Hu Jie. Three dimensional spectral-spatial feature 
extraction and small sample classification technology for Hyperspectral 
remote sensing image [D]. Shenzhen: Shenzhen University, 2017. ) 

[16] Grupo De Inteligencia Computacional. Hyperspectral remote sensing scenes 

[EB/OL]. (2014-04-07) [2018-05-16]. http://alweb. ehu. es/ccwintco/index. 


php?title=Hyperspectral Remote Sensing_Scenes. 


